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5 : ABSTRACT 



We calculate temporal correlators and their spectral functions with meson quan- 
tum numbers in the deconfined phase of QCD using the hard thermal loop (HTL) 

approximation for the quark propagator. Although this approach does not result 
• 1—1 , 

^ ■ in a complete next-to-leading order perturbative calculation it takes into account 

H ; 

5t , important medium effects such as thermal quark masses and Landau damping in 

the quark-gluon plasma. We show that both effects lead to competing modifications 
of the free mesonic correlation functions. We find that correlators in scalar channels 
are only moderately influenced by the HTL medium effects, while the HTL-vertex 
corrections lead to divergent vector correlators. 
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1 Introduction 



While the lattice calculations of hadron properties in the vacuum have reached quite 
satisfactory precision, little is known from such first principle calculations about 
basic hadronic parameters in a thermal medium, e.g. masses and widths at finite 
temperature. Lattice calculations of such quantities at zero temperature generally 
proceed through the calculation of correlation functions in Euclidean time. This 
approach naturally carries over to the calculation of spatial correlation functions 
at finite temperature. Such spatial correlation functions indeed show evidence for 
sudden changes of in medium hadron properties above Tc 0. However, they provide 
only indirect evidence for modifications of, e.g. hadron masses and their widths. 
The appropriate approach here would be a detailed analysis of temporal correlation 
functions 0, |^, which at finite temperature are restricted to the Euclidean time 
interval [0, 1/T]. The interesting information on hadronic states is then encoded in 
the spectral functions for these correlators [Q. 

Currently available results from lattice calculations show significant changes in 
the behaviour of temporal correlation functions in the high temperature plasma 
phase of QCD H, H, ^, H. However, at least close to Tc the correlation functions 
clearly deviate from those of freely propagating quarks. It thus is important to 
understand in how far the temporal correlation functions carry information about 
the existence or non-existence of bound states or resonances in the plasma phase. 
Various calculations within the framework of low energy effective models also suggest 
strong modifications of hadron properties and consequently also of the spectral 
functions [§]. However, it is difficult in such model calculations to deal with the quark 
substructure of hadrons, which will become important at high temperature where 
one expects to find indications for the propagation of almost free, massless quarks. 
Eventually it is the hope, that spectral methods 0, which successfully have been 
applied to hadron spectral functions at zero temperature [|lOl, can also be applied 



at finite temperature. In particular, in the high temperature limit, well above the 
QCD phase transition temperature, it then might be appropriate to compare lattice 
calculations for temporal hadron correlators also with perturbative calculations. At 
least to some extend non-perturbative information can also be incorporated in such 
an analytic calculation by using the hard thermal loop (HTL) resummation scheme 
Tl[] . The recent successes in reproducing the QCD equation of state calculated 



on the lattice [12, O] with HTL-resummed perturbative calculations for T > 2Tc 



T^, |T5l suggest that this may be a reasonable starting point also for the description 



of other properties of the high temperature plasma phase [0. Indeed, the spectral 
functions, which one will extract from an analysis of temporal correlators, are closely 
related to quark- ant iquark annihilation processes in the quark-gluon plasma. For the 
vector channel this is linked to the dilepton production at high temperature, which 



has been studied in the HTL-approximation The temperature dependence of 
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the pseudo-scalar correlator is related to the chiral condensate. Thermal effects 
on this as well as the pseudo-scalar masses and dispersion relations inffuence the 
appearance or suppression of a disoriented chiral condensate which might lead to 



observable effects in relativistic heavy ion collisions [IS 



We will analyze here the structure of temporal correlation functions within the 
context of HTL-resummed perturbation theory. In the infinite temperature limit 
the free field behaviour is expected to give the dominant contribution to the spec- 
tral functions for energies ~ T |T0, 120]. Here we will use the HTL-resummed 



quark propagator ^ for calculating the temporal correlators and their spectral 
functions. In this way important in-medium properties of quarks in a quark-gluon 
plasma are taken into account]^. As we will see, these medium effects, will lead 
to competing effects in thermal correlation functions. On the one hand the effec- 
tive quark propagator takes into account the generation of thermal quark masses, 
rriT ~ g{T)T. This cuts off the low frequency part in the spectral functions and thus 
will lead to a steepening of thermal correlation functions. On the other hand it also 
contains the contributions from plasmino modes as well as interactions of quarks 
and antiquarks with gluons in the thermal heat bath (Landau damping). This en- 
hances the contribution of soft modes with u ~ g(T)T which, in fact, will dominate 
the structure of spectral functions at low energies even at rather high temperature. 
These contributions from soft modes will lead to a flattening of thermal correlation 
functions. We will discuss the interplay between both features of HTL-resummed 
correlation functions in this paper. 

In the next section we will present the framework for the calculation of thermal 
meson correlation functions in the HTL-approximation and give results for the scalar 
and vector spectral function. In Section 3 we compare the resulting thermal meson 
correlation functions with the leading order perturbative (free) correlators. Finally 
we give our conclusions in Section 4. 



2 Thermal Meson Correlation Functions 



2.1 Definitions 

We want to analyze the behaviour of meson correlation functions in the high temper- 
ature limit. They are constructed from meson currents Jm{j,x) = q{T,x)TMq{'T,x), 
where Tm is an appropriate combination of 7-matrices that fixes the quantum num- 

^It should be noted that we are not aiming at a calculation in the strictly perturbative sense, 
which holds only in the weak coupling limit, g ^ I. Rather we want to consider QGP medium 
effects using the HTL quark propagator for comparing with lattice results at realistic values of the 
coupling constant. 
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bers of a meson channel; i.e., Tm = 1, 75, 7/^, 7/^75 for scalar, pseudo-scalar, vector 
and pseudo-vector channels, respectively. The thermal two-point functions in coor- 
dinate space, Gm{t,x), are defined as 



GM{r,x) = (JAf(r,f)jl,(0,0)) 

= T Y: / 7^e-^^^"^""^^XMK,p) , (2.1) 

where r G [0, 1/T], and the Fourier transformed correlation function Xui'^niP) is 
given at the discrete Matsubara modes, ujn = 2n7cT. The imaginary part of the 
momentum space correlator gives the spectral function (Jm{'^i'P)i 

XM{uJn,p} = - f duj . ^^^^'^ =^ (rM{i^,P) = -Im Xm{^,P) ■ (2.2) 

J -00 lUOn — 00 + te TC 



Using eqs. ^3] and |2.2| we obtain the spectral representation of the thermal cor- 



relation functions in coordinate space at fixed momentum {(3 = 1/T), 
n t ^ r A ( ^ cosh(a;(r-/3/2)) 

GM{r,p)= dujaM{uj,p) . ^ • 2.3 

JO smh(a;p/2) 

For a perturbative analysis of these correlation functions at high temperature it 
is convenient to introduce dimensionless variables, uj = u/T, p = p/T, f = tT, 
the reduced spectral function a{Lj,p) = a{uj,p)/T'^, and the reduced correlator 

GM{fJ)=GM{T,p}/T^ 



2.2 Free Meson Spectral Functions 

The starting point for a calculation of the meson spectral functions and the me- 



son correlation functions is the momentum space representation of the latter pO |. 
To leading order perturbation theory one has to evaluate the self-energy diagram 
shown in Fig. la, where the internal quark lines represent a bare quark propagator 
S'i?(/co,fc) which can be expressed in terms of its spectral function pp{u,k,m) and 
is conveniently written as 

Sp{ko, k) = -(70 ko-rk+m) / dre'^^^ / duppiu, k,m) [l-npiuj)] e"'^^, (2.4) 

Jo J -00 

where ko = {2n + l)z7rT, npiuj) = 1/(1 + exp(u;/T)) and 
1 

PF{u;,k,m) = — {5{uj - Uk) + 5{uj + ujk)) , (2.5) 
2uj 





Figure 1: The self-energy diagrams for free quarks (a) and in the HTL approximation 
(b). 



with iUj. = y k"^ + m? 



The thermal meson spectral functions are then obtained from eq. |2.2| , that is 
from the imaginary part of the correlation functions in momentum space, 



TMSpiko, k)T\jSy{uj - ko,p- k) 



(2.6) 



In the case of free fermions this is easily evaluated. In the limit of vanishing external 
momentum one finds for the spectral functions, 



47r2 



e(c<; - 2m) tanh(Lj/4r) 



, 2m\2 



(2.7) 



where different quantum number channels are characterized by the pair of param- 
eters {aMjhu)- For the scalar (s), pseudo-scalar (ps), vector (v) and pseudo-vector 
(pv) channels they are given by (-1,1), (1,0), (2,1) and (-2,3), respect ively[|. In 
the massless limit the spectral functions are chirally symmetric, \(Jps\ = \(ys\ and 
|cr„|. In this case the remaining integral in eq. p.3| can be done analytically 



pv I 



and one obtains for example in the pseudo-scalar case ||2^ 



G'p.(f,p = 0) = 27riV, (1-2 



^ 1 + cos2(27rr) 
^' sin3(27rf) 



+ 4iV, 



cos(27rr) 
sin^(27rf) 



f2.8) 



2.3 HTL- Approximation for Meson Spectral Functions 

Now we want to go beyond the free quark approximation and consider in-medium 
quark propagators. A consistent way in the weak coupling limit (^f <^ 1) is the 

■^In the vector and pseudo-vector cases we denote by crj\/ the trace over the Lorentz indices of 
111/ 
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use of a HTL-resummed quark propagator (and quark-meson vertex) if the quark 
momentum is soft, i.e. of order gT. These corrections to the spectral function are 
of the same order as the free spectral function for small energy u. 



Using the HTL resummation technique important medium effects of the quark- 
gluon plasma such as effective quark masses and Landau damping are taking into 
account. The HTL-resummed fermion propagator is obtained from eq. ^.4| by replac- 
ing the free spectral function pp with the HTL-resummed spectral function which 
for massless quarks is given by []T7|, pT| 



PHTL(fco, k) = ^p+{ko, A;)(7o - 2 A; ■ 7) + ^P-(fco, k){'yo + ik- 7) 



(2.9) 



with k = k/k, k = \k\, and 



P±{kQ,k) 
P±{ko, k) 



2m^ 



m: 



T 



[S{ko - u;±) + 6{ko + u;^)] + P±{ko, k)Q{e - k^ 



±kn - k 



k{-ko ±k) + ±1 - ±^ ln|±g 



(2.10) 



TT 2 ±ko-k 
2 '"-T k 



Here u}±{k) denote the two dispersion relations of quarks in a thermal medium 
T7| , |21[] and uit = g(T)T/\/6 is the thermal quark mass. We note that in addition 



to the appearance of two branches in the thermal quark dispersion relation the 
HTL-resummed fermion propagator also receives a cut-contribution below the light- 
cone (fcg < k'^), which results from interactions of the valence quarks with gluons 
in the thermal medium (Landau damping). Furthermore, an explicit temperature 
dependence only enters through mT{T). Also the HTL-resummed quark spectral 
function can thus be written in terms of dimensionless, rescaled variables, e.g. uj = 
uj/T etc. and the reduced meson spectral functions 

~HTL ^ ^HTL 

/T^ will depend 

on temperature only through •ffij' = g{T)/\/Q. It also should be noted that the 
HTL resummed quark propagator is chiral symmetric in spite of the appearance of 



an effective quark mass ||2T[. In the following we thus will ignore the parity of the 



meson states and will generically talk about scalar and vector channels onlyQ. 



Inserting eq. |2]^ for m = together with eq. ^]9| and eq. p.lO| into eq. |2]^ we 



can determine the spectral functions for mesons in the HTL-approximation. In the 
vector channel this also requires additional modifications of the vertex functions Fm, 
i.e., the use of a HTL quark-meson vertex, as discussed in detail in fl^. In the case 



of the scalar and pseudo-scalar spectral function the vertices are given by the bare 



■^We will, however, show results for pseudo-scalar and vector spectral functions and correlators 
which in our notation are strictly positive. 
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vertices = 1 and Tps = 75, since contributions from a HTL resummation to the 
vertices are suppressed in this case, i.e. lead to higher order corrections, as discussed 
for the scalar case in the Yukawa theory ||2^ and scalar QED In the case of 



QCD the absence of HTL corrections for scalar as well as pseudo-scalar vertices has 
been shown in the context of kinetic theory p4|. If, on the other hand, the bare 



vertex is proportional to 7^ as in the vector meson case, the HTL vertex cannot be 
neglected, since in a gauge theory it is related to the HTL fermion propagator by 
Ward identities 



The pseudo-scalar spectral function can then be written as 
(Tps{uj,p} = 2Nc{e'^''^ - 1) / -— - / dxdx' nF{x)nFix')5{uj - x - x') 

J (zTT)'^ J-oo 



(1 - q- k)[p+{x,k)p+{x',q) + p-{x,k)p_{x,q)] 
+ (1 + q- k)[p+{x, k)p^{x', q) + p_(x, k)p+{x, (2.11) 



where q = p — k. The corresponding relation for the vector spectral function, which 
also includes HTL-vertex contributions, is related to the dilepton production rate 
calculated in Ref. |]T^ in the HTL approximation. 



a.M= 0) = ^ (e-/- - 1) -^{p = 0) . (2.12) 

Here a is the electromagnetic fine structure constant. 

As the thermal meson correlation functions are constructed from two quark prop- 
agators, they will receive pole-pole, pole-cut and cut-cut contributions, i.e., the 
mesonic spectral functions for p = are generically given by 

a^^^iu) = (tPp(cj) + ctP=(cj) + a^^{uj) . (2.13) 

Explicit expressions for the three different contributions to the pseudo-scalar spectral 
function, (J^^^, are given in the Appendix. Similar results for the vector spectral 
function a^^^ have been derived in |T^ where also a detailed discussion of the 



physical processes related to the pole-pole, pole-cut and cut-cut contributions is 
given. In particular, there are characteristic peaks that show up in the pole-pole 
contribution (van Hove singularities). They are caused by a diverging density of 
states which is inversely proportional to the derivative of the dispersion relations, 
uj'^{k), appearing in eq. A.l. Owing to the minimum in the plasmino branch^ these 

®In Ref.[^ it has been argued that the full in- medium quark propagator leads in general to 
two branches in the dispersion relation, of which one exhibits a minimum. 
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Figure 2: The pole-pole, pole-cut and cut-cut contributions to the pseudo-scalar 
(a) and vector (b) spectral function for tjit = 1- The crosses show the free meson 
spectral function. 



derivatives vanish at = 0A7mT and 1.856mT. Apart from values close to the van 
Hove singularities in 0"^^ one finds that the cut contributions dominate the spectral 
function for small values of a), e.g. for u < g{T). 

In Fig. 2 we show the pole [a^^) and cut (5"'"^, a^'^) contributions to the scalar 
(Fig. 2a) and vector (Fig. 2b) spectral functions for the case rfiT = 1, extrapolating 
the HTL results, obtained in the weak coupling limit, to g = \/6. Of course it is a 
priori not clear whether one can apply the HTL resummation method to such large 
couplings, although the extrapolation to realistic values of the coupling constant 
has been used for various observables of the QGP However, in this way we 



can study the influence of medium effects at least qualitatively and as will become 
clear, the small corrections found by us do justify our choice for the couplings. As 
can be seen in Fig. 2 the pole and cut contributions influence the spectral function 
in different ways. The former gives the dominant contribution for large cD. The de- 
viations of cr^"""^ from the free spectral function in this energy regime as well as the 
threshold for ci) ~ 2 rfiT is due to the presence of a non-vanishing thermal mass in 
the quark dispersion relation and reflects the almost free propagation of two quarks 
in the plasma. Additional interactions of these quarks with the thermal medium 
(Landau damping) are represented by the cut contributions. These lead to an en- 
hancement over the free spectral functions for small values of u as discussed above. 
Furthermore, we note that the pole-pole and pole-cut contributions to the spectral 
functions are similar in the scalar and vector channels. The cut-cut contribution, 
however, behaves differently at small energies. While it vanishes for small Cj in the 
scalar channel it diverges linearly in the vector channel. This can be traced back to 
the structure of the effective HTL- vertex, which contains a coUinear singularity 



As a consequence of this singularity infinitely many higher order diagrams in the 
HTL expansion contribute to the same order in the coupling constant |2^. This, 
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Figure 3: The thermal pseudo-scalar meson correlation function in the HTL ap- 
proximation for friT = 1 (left) and rfiT = 2 (right). The curves shown the complete 
thermal correlator (middle line), the correlator constructed from (jPP only (lower 
line) and the free thermal correlator (upper line). 

of course, indicates that the low frequency part of the vector spectral functions is 
inherently non-perturbative. 

3 HTL- Approximation for Thermal Meson Cor- 
relators 

3.1 Thermal Pseudo-Scalar Meson Correlation Function 

Now we will compute the temporal correlators from the spectral functions, derived 
in the last section, using (|2.3| ). Although the temporal correlators are dominated 
by hard energies ~ T, we use the spectral functions, calculated within the HTL 
resummation scheme, in order to take into account medium effects from the quark 
propagator at least semi-empirically. 

The competing influence of pole and cut contributions to the HTL-resummed 
spectral functions carries over to the behaviour of thermal meson correlation func- 
tions. The appearance of a non-vanishing thermal quark mass tends to lead to a 
more rapid decrease of the correlator in Euclidean time than this is the case for the 
free massless correlator. The enhancement of the low energy part which is due to 
the cut contributions, on the other hand, will counteract this trend. This is evident 
from the behaviour of thermal meson correlation functions in the scalar channel 
which is shown in Fig. 3 for rfiT = 1 and 2. As expected the correlator constructed 
from the pole-pole contribution alone is steeper than the free correlator and, more- 
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Figure 4: The ratio of the HTL-resummed and free thermal pseudo-scalar correlation 
function versus Euclidean time r in units of the temperature. Shown are results for 
thermal quark masses mr = 0.5 (top), 1.0 (middle) and 2.0 (bottom). 



over, is strongly dependent on ttit- The cut contributions, however, enhance the low 
energy contributions in the spectral function and thus flattens the correlator again. 
Somewhat surprisingly for rfiT — 1 this seems to compensate almost completely the 
deviations from the free correlator introduced by the pole contributions. The differ- 
ence between the free and HTL-resummed correlators is largest for tT ~ 1/2 where 
the contribution from the low energy regime in the spectral function is largest. For 
tT 0, and tT ^ 1, on the other hand, the free and HTL-resummed correlators 
approach each other as lim^^^oo o"^"'"'"(^)/'7^'^'''^(^) = 1- These features are amplified 
in the ratio G^^^ {r) / Gp^'^ {r) which is shown in Fig. 4. 



3.2 Thermal Vector Meson Correlation Function 



As already discussed in section 2.1 the calculation of the vector correlators within 
the HTL method requires the use of effective quark-meson vertices as shown in 
Fig. lb. This does lead to a hnear divergence of the spectral function in the vector 
channel at low frequencies, which in turn renders the temporal correlator infrared 
divergent. In fact, although the scalar correlation functions are infrared finite, it is 
to be expected, that also in this case the low frequency part of the HTL-resummed 
spectral functions will be modified significantly from contributions of higher order 
diagrams. It thus seems to be reasonable to consider modified correlation functions, 
which are less sensitive to details of the low frequency part of the spectral functions. 
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Figure 5: The subtracted thermal pseudo-scalar (left) and vector (right) meson 
correlation functions in the HTL approximation for rfiT = 2. The curves shown the 
complete thermal correlator (middle line), the correlator constructed from a^f only 
(lower line) and the free thermal correlator (upper line). 

We therefore define the subtracted correlators 

AGm{t) = Gm{t) - Gm{P/2) . (3.1) 

In the subtracted correlation functions the infrared divergences are eliminated. They 
arc well-defined in the scalar as well as in the vector channels. In Fig. 5 wc compare 
the HTL-resummed subtracted correlation functions with corresponding results for 
the free case. This shows that after ehmination of the infrared divergent parts the 
structure of the pole and cut contributions is similar in scalar and vector channels. 
The vector correlator seems to be even closer to the leading order perturbative (free) 
correlator than the scalar correlation function. 



4 Conclusions 

We have calculated thermal meson correlation functions and their spectral functions 
in the HTL-approximation; i.e., we have used the HTL-resummed quark propagator 
(and HTL quark-meson vertex in the vector case) for considering important in- 
medium properties of quarks in a QGP. 

Wc have analyzed the infiuence of the various contributions to the HTL-resummed 
scalar and vector meson spectral functions on the structure of the thermal corre- 
lators. We generally find that the correlators in HTL-approximation - after sub- 
tracting the infrared singularity in the vector correlator - are quite similar to those 
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calculated in leading order perturbation theory which correspond to free correlation 
functions. Non-perturbative features of the HTL-resummed quark propagators such 
as the generation of a thermal quark mass and Landau damping are clearly visible 
in the meson spectral functions. However, they lead to competing effects in the 
correlation functions and to a large extent compensate each other. 

The main difference between the scalar and the vector channel using the HTL 
approximation is the different behaviour of the cut-cut contribution to the spectral 
functions at small energies. The vector spectral function diverges in the infrared 
limit leading to a singular expression for the vector correlation function. This feature 
is in agreement with the observation, that the dilepton production rate, which is 
closely related to the vector spectral function, cannot be computed within the HTL 
improved perturbation scheme for small invariant masses M ~ g^T. 

The existing lattice calculations of thermal meson correlation functions show that 
the correlators deviate from the free field result significantly for temperatures T<2Tc. 
However, also at larger temperatures it seems that the scalar correlator only slowly 
approaches the free correlation function and still differs in shape from the vector 
correlator. This suggests that HTL-resummed perturbation theory, which gave a 
satisfactory description of bulk thermodynamics above 2Tc, will not be appropriate 
for a quantitative analysis of thermal meson correlation functions at least in the 
pseudo-scalar case. In other words, the HTL medium effects (thermal quark masses. 
Landau damping) are not sufficient to explain the deviations of the pseudo-scalar 
correlator from the free one as observed in lattice calculations. Therefore additional 
non-perturbative effects, maybe related to chiral symmetry restoration, appear to 
be important in the pseudo-scalar channel. 
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Appendix 

We will give here explicit expressions for the pole-pole, pole-cut and cut-cut contri- 
butions to the pseudo-scalar spectral function. The pole-pole contribution (p — 0) 
is given by 




27r2 



rrirp 
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i=l 



W^{k\) - Uj'_{k\)\ 



i=l 



2\uj'_{kl)\ 



(Al) 



Here uj±{k) denote the quark dispersion relations for the ordinary quark (+) and 
the plasmino (-) branch ki is the solution of u — 2u^{ki) = 0, ^2 ci-nd /C3 are the 
solutions of u; — (^+(^2) + u!_{k2) = and u — 2uj^{k\) = 0, respectively. Note that 
for small momenta the last two equations can each have two solutions. Furthermore, 
uj'_^{k) = {duj±{x) / dx)\x=k- For the pole-cut contribution we find 



2Nr, 



TT 



2 "^T 



2(e-/r_i) r^k k^ 

Jo 



+ 9(A;^ — (u; — uj^)'^)np{uj — uj_)np{uj_)P_{uj — uj_, k){uj'^ — k"^) 

+ e(A;^ -{00 + uj_ f)np{uj + uj_)[l - np{uj_)]l3+{uj + uj_, k){uj't - k'^) 

+ e{P -{uj + uj+f)np{uj + cj+)[l - np{uj+)]p+{uj + iu+, k){tul - /c^) 

{A.2) 



Finally we obtain for the cut-cut contribution 



2N, 



^;^(e"/^ - 1) / dkP dx np(x)np(x - uj)Q(e - (x - uf) ■ 

TT^ JO J-k 

P+{x, k)l3^{uj — x,k) + P-{x, k)P_{uj — X, k) 

(A3) 



The resulting contributions to the thermal correlator are then given by eq. p73 
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